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Abstract 

We construct a hadronic equation of state consistent with chemical freezeout and discuss how 
such an equation of state modifies the radial and elliptic flow in a hydrodynamic + hadronic 
cascade model of relativistic heavy ion collisions at the SPS. Incorporating chemical freezeout 
does not change the relation between pressure and energy density. However, it does change the 
relation between temperature and energy density. Consequently, when the hydrodynamic solution 
and freezeout are expressed in terms of energy density, chemical freezeout does not modify the 
hydrodynamic radial and elliptic flow velocities studied previously. Finally, we examine chemical 
freezeout within the hadronic cascade (RQMD). Once chemical freezeout is incorporated into the 
hydrodynamics, the final spectra and fireball lifetimes are insensitive to the temperature at which 
the switch from hydrodynamics to cascade is made. Closer inspection indicates that the pion 
spectrum in chemically frozen hydrodynamics is significantly cooler than in the hydro+cascade 
model. This difference is reflected in V2{pt)- We extract the freezeout hadron density in RQMD and 
interpret it in thermal terms; the freezeout hadron density corresponds to a freezeout temperature 
of Tf « 100 MeV and /% ps 80 MeV. 
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I. INTRODUCTION 



Currently, by colliding heavy ions at the Super Proton Synchrotron (SPS) and the Rel- 
ativists Heavy Ion Collider (RHIC), experimentalists [|l], |2| have endeavored to create a 
deconfined state of quarks and gluons - the Quark Gluon Plasma (QGP). Some degree of 
thermalization in the heavy ion reaction is a prerequisite for QGP formation. It is an ex- 
perimental fact that the ratios of the twenty or so hadron species produced in the heavy ion 
collision are close to the thermal ratios expected of an ideal gas of hadrons at a temperature 
of T ?a 160 — 180 MeV ||. It is exciting that this temperature is close to the transition tem- 
perature to the QGP, T c pa 160 MeV. However, after accounting for finite size corrections, 
the hadron ratios in pp and e + e~ reactions are reproduced by the same thermal models 
used to describe heavy ion data. It seems that statistical mechanics provides a universal 
description of hadronization. Unlike pp and e + e~ reactions, the ~ 5000 hadrons produced in 
the heavy ion collision rescatter after hadronization. Ideally, this ensemble of hadrons may 
be considered a hot hadronic gas; hydrodynamics then describes the subsequent evolution. 
With an understanding of the final hadronic expansion, the global properties of the collision 
can be quantified. 

We first review the notion of chemical freezeout in the hadron gas [f|, |5|, ^ [Tj. In the 
hadron gas, the principal hadronic reactions, e.g. Tin — ► p — ► Tin, ttK — > K* — > 7r.fr, 
ttN — > A — > ttN, do not change the yields of pions, kaons, and nucleons. Nevertheless, these 
reactions provide a mechanism for thermalization during the hadronic evolution. Thermal 
equilibration times at a temperature of T ~ 160 MeV are typically ~ 2fm/c ||. In the 
strong interactions of light hadrons there are only three conserved currents: baryon number 
(</#), strangeness and isospin (Jf). Other hadronic reactions (for instance irN — > 

iV*(1530) — > An — > Ntttt) do change the total yield of pions, kaons and nucleons. Because 
of these reactions the system approaches chemical equilibrium. The yield of pions, kaons, and 
nucleons changes via such reactions until the Gibbs free energy reaches a minimum. However, 
the time scale of chemical equilibration is much longer than that of thermal equilibration; 
chemical equilibration times are typically r c h ~ 200fm/c ||, rO | . Therefore in the hadron 
gas, there are two disparate time scales, r t h and r c h- 

In heavy ion collisions, the lifetime of the hadronic stage is approximately ~ lOfm/c, 
which is very short compared to r c h , but longer than Tth- Therefore, on the time scale 
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of the collision, although the chemical composition is fixed at the time of hadronization, 
the system continues to evolve kinetically for some time until the particles breakup. The 
stages of the collision have been described with the following picture. First, there is chemical 
freezeout (hadronization). Then the hadrons evolve as a hadronic fluid until thermal freeze- 
out (breakup). In this picture, the particles develop chemical potentials during the hadronic 
evolution since the total number of particles is fixed and the temperature decreases. 

There is some evidence for this chemical/thermal freezeout picture. First, hadronic cas- 
cades indicate that following hadronization pions and nucleons rescatter for ~ 10 fm/c. 
Consequently, the pion (pt) decreases while the nucleon (pr) increases. Thus, microscopic 



calculations indicate that pions cool and increase the transverse flow of the nucleons [ I I 



Second, Rapp and Shuryak have argued that a pion chemical potential can explain the 
anti-baryon yield at the SPS, since p and p can be produced via the forward and backward 
reactions of pp <-> 5ir |12[ . Finally, the most compelling experimental evidence for the pion 



chemical potential comes from a combined analysis of spectra and HBT correlations [13 



This study indicates that the hadron density at freezeout is substantially lower than the 
density of a hadron gas at T « 160 MeV. Furthermore, the pion phase space density, which 
can be extracted from from two pion correlations, is overpopulated. An overpopulated phase 
space is expected from a Bose gas with a chemical potential. 

Hydrodynamics has been used extensively to model heavy ion reactions. However, these 
hydrodynamic calculations assume thermal and chemical equilibrium. While thermal equi- 
libration is at least plausible, chemical equilibration is certainly impossible. Typically a hy- 
drodynamic simulation is run until a universal freezeout temperature Tf « 120 — 140 MeV, 
which is adjusted to match the pion and nucleon spectra. Subsequently, even if the spectral 
shape is correct, the yields of p,A,K for example, are typically wrong by factors of two. 
To account for this discrepancy, a comparison to the data is then made by rescaling the 
particle yields (by hand) to their value at T « 160 MeV. Energy-momentum and number 
conservation are violated in this inconsistent procedure. 

One approach to the problem of chemical freezeout is to stop the hydrodynamic evolution 
at T c and then continue the evolution with a hadronic cascade JIT], fnj. This approach 
incorporates chemical freezeout as chemical equilibration times are encoded into the hadronic 
cross sections in the cascade. Although hydro+cascade provides a comprehensive description 
of the heavy ion data and has been used extensively, the final spectra are sensitive to exactly 
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where the switch to the cascade is made JTT|]. Ideally, the switching temperature, T SW i tc h, 
should be varied, with final results independent of this artificial parameter. However, the 
cascade conserves n, K, A . . . number in the dominant reactions; therefore to achieve a kind of 
dual description chemical freezeout must be incorporated into the hydrodynamic evolution. 

This incorporation may be achieved by including additional conservation laws for T < T c . 
For each conserved species it, K, A • • •, we have 



0. 



(1) 



The Equation of State(EOS) is modified and now depends on n n , Uk, ■ • '> m addition 
to e and ub- Furthermore, the relationship between energy density and temperature is 
dramatically different. In this work, we consistently incorporate chemical freezeout into 



the hydrodynamic evolution of the hydro+cascade model, H2H fL4[ . Similar recent efforts 
have been presented J15, 16 1. H2H [14] and other hydrodynamic models [17, 18] have been 
compared extensively to available data both at the SPS and RHIC. The purpose of this work 
is to illustrate how chemical freezeout changes the results of these works. 

Briefly, in the H2H model the initial stages of the collision are modeled with (2+1) ideal 
fluid dynamics, assuming Bjorken scaling in the longitudinal direction [TS , BIJ. The initial 
entropy is distributed in the transverse plane according to a Glauber model at a time, 
t — 1 fm/c. At particular temperature, T switc h (which is less than the critical temperature 
T c = 165 MeV), the fluid is converted into hadrons via the Cooper- Frye formula [21]. The 
Cooper-Frye formula is appended with a theta function rejecting backward moving hadrons. 
Initial conditions for the hydrodynamic evolution have been chosen to match the model 
charged particle multiplicity and net baryon number to the experimental values at the SPS 
(PbPb with v/i = 17 GeV A) and RHIC (AuAu with y/s = 13Q GeV A) |4|. 

The calculation presented here are for the SPS (PbPb with y/s = 17 GeV A) at an 
impact parameter of b=6 fm. (A non-central impact parameter is taken in order to study 
elliptic flow. Final flow velocities and freezeout conditions do not change rapidly with impact 
parameter |TJ|, j22f.) The initial conditions (see Ref. |T^| for more detail) are characterized by 
an entropy per baryon ratio, s/ub = 42, and an average initial energy density, < e >| To = 
5.4 GeV/fm 3 . The EOS has an 0.8 GeV/fm 3 latent heat and is identical (above T c ) to 
the equation of state LH8 which was used previously |14|]. As this work concentrates on 
the freezeout stages of the collision, the differences between the presented results and the 
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analogous calculations at RHIC are small | |23|| . The slightly larger baryon density at the 
SPS does not significantly alter the freezeout dynamics which is already meson dominated 
at the SPS 0, |TJ 

Sect. [IT] discusses the EOS in which all the hadronic species are conserved. The hydrody- 
namic equations are then solved with and without the additional conservation laws and the 
solutions contrasted. The principal result is that the particle yields of a hydrodynamic cal- 
culation can be consistently modified (increased), provided the freezeout temperature is also 
modified (decreased), to keep the energy density the same. Returning to the hydro + cascade 
approach, in Sect, [ill], we vary the switching temperature to the cascade, T switc h- For 
115 MeV < T switc h < 160 MeV, the model spectra are insensitive to variations. Elliptic 
flow remains somewhat sensitive to the switching temperature, although the sensitivity is 
reduced when chemical freezeout is incorporated. 

With these basic results, the freezeout conditions of the cascade RQMD are analyzed in 



Sect. UJ. At freezeout, the density of pions in the cascade is equal to the density of an ideal 
hadron gas with Tf « 100 MeV and /x w ~ 80 MeV. This freezeout density is universal and 
is independent of impact parameter, collision energy and the switching temperature, T SW i tc h 
fTT| , |Tj| . In evaluating the extent to which the cascade reproduces the hydrodynamics in 
this low temperature region, we find that the hydrodynamics cools much more quickly than 
the cascade. This rapid cooling impacts both spectra and v 2 (pt)- 



II. IMPLEMENTING CHEMICAL FREEZEOUT IN THE EOS 

In this section, we construct an EOS which incorporates chemical freezeout. (For similar 
constructions see @, |K| [OJ j24|.) Above T c , the EOS is identical to LH8 which was described 
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14| . Only the hadronic portion of the EOS is modified. 
For the purposes of this work, we consider expanding hadronic gas with only the lowest 
hadron multiplets. Specifically, we consider only the 0~ and 1~ meson octets, rj' and 0, the 
i + baryon and anti-baryon octets and the | + baryon and anti-baryon decuplets. The strong 
resonant reactions in this ensemble excite mesons and baryons from the lower hadronic mul- 
tiplet to the higher multiplet and are given (up to isospin and baryon/anti-baryon symmetry) 
by 



7T7T 



P — > 7T7T (2) 



7lK + -> K + * -> 

vriV -> A -»■ vriV 
txY -> E* -> ttF 



where Y denotes the hyperons A, E. These are certainly the most important reactions for 
the late stage of the heavy ion collision. If only these reactions and elastic collisions are 
included, then the system has 16 conserved currents: baryon number (Jg), strangeness( Jg), 
isospin (Jf) and 13 other conserved numbers, J#., where Hi runs over the hadron species, 

n,N,Y, Y, S, S, K~,r), ?/, uj, 0, fi, 0. (3) 

Nucleon number and number are not included in the list of hadrons, but their conser- 
vation follows from the conserved hadron currents already specified and from baryon and 
strangeness conservation. Perhaps some of the stranger species (Q, (f>) should not be included 
in the list of hadrons since they decouple early in the hadronic evolution. In addition, u 
does not neatly fit into the thermal/chemical freezeout picture, as its lifetime is comparable 
to the collision lifetime. Nevertheless, these particles are included for theoretical consistency 
and we have found that they do not affect the bulk properties. 

With the assumption of equilibrium the hydrodynamic equations become 

d^ u = (4) 
d^mn") = 0, (5) 

where T^ v = (e + ■p)u il u v — pg^". This set of equations, together with the assumption of 
smooth flow, imply entropy conservation, 

d^sm*) = 0. (6) 

An EOS - a relation between the pressure and the energy density and number densities 
associated with the 16 conserved currents - is required to complete the system. 

The hadronic EOS is taken as a sum of independent ideal gases over the hadrons con- 
sidered. This approximation is based upon the fact that the thermodynamics of pions 
interacting via irn — > p — > nn is nearly equivalent to the thermodynamics of an ideal gas of 
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7r-s and p-s |[25|| . In all the ideal gas formulas, the 16 chemical potentials all appear in the 
combination, 

jjt, B B + fi s S + p/J + VmHi, (7) 

i 

where Hi runs over the hadron list given in Eq. [| For instance, 7r~ has / = —1/2, if,,- = 1, 
and all other quantum numbers zero. On the other hand, K~* has B = 0, S = 1, / = —1/2, 
iJ^ = 1, H K - = 1 and all other quantum numbers zero. The EOS is a relationship between 
the pressure p, and the energy density e, and the number densities {rii}, and is therefore a 
function of 1+16=17 variables. 

In the heavy ion collision, only a small fraction of the total 17 variable phase space is 
sampled. Indeed, if at some time n,/s is constant everywhere in space, then Eq. |5| and 
Eq. |] imply that rii/s is a constant for all later times [26|, jZ7fl . (This follows directly since 



9 M log(s) = <9 M log(rij) = —d^ .) Thus with the assumption of entropy conservation, it is 
sufficient to know the EOS along the trajectory where rii/s is constant. We now specify 
rii/s and construct the hadronic EOS along this trajectory. This unique trajectory in the 
17 variable phase space is called the adiabatic path below. 

The ratio s/ub = 42 is fixed by the initial conditions and is chosen to reproduce the 
experimental proton to pion ratio ||14j| . The net strangeness is zero, n s = 0. For simplicity, 
the net isospin is taken to be zero, rij = 0. As discussed in the introduction, at hadronization 
the system is born into chemical equilibrium. Therefore at T c , tih/s is taken to be its 
chemically equilibrated value. Now, with n^/s a specified constant on the T = T c hyper- 
surface, as the system expands and cools, rii/s remains constant. The particles develop 
chemical potentials to ensure this constancy. At a computational level, the procedure for 
constructing the EOS along the adiabatic trajectory is the following. First, at T = T c adjust 
liBi^s-i an d fJ-i, until ns = nj = and s/ns = 42. Then, calculate n^/s at T c . Then in small 
increments lower the temperature, adjust all the chemical potentials to leave n,/s constant, 
and tabulate all the thermodynamic quantities (i.e, pressure, entropy, energy density) along 
the way. 

Consider the results of this procedure. First, a few of the chemical potentials are shown 
as a function of the switching temperature by the thick lines in Fig. |l|(a). We can derive 
an approximate formula for the chemical potentials as a function of T, for a collection of 
non-relativistic ideal gases. For non-degenerate and non-relativistic ideal gases the partial 




FIG. 1: Chemical potentials as a function of temperature at the SPS (s/ns=42) for (a) ir, K and 



A and (b) it only. Eq. 11 gives an approximate formula for these chemical potentials. 



pressure pi, and partial energy density of the i-th species are given by, 

Pi = n{T (8) 

&i = riiiTii + (Const)njT. (9) 

Since T s = J2i^i+Pi~ Vi^i , 

l = E=i(=^ + Const). (10) 



s V T 

rii/s is constant below for T < T c . To keep all the terms in parenthesis constant below T c , 
we require 

Mi = ™* (l - 5^) + ^FT, (11) 

where fi c is the chemical potential at the critical temperature. The approximation is shown 
by the dashed lines and works well for all particles except pions. Being the most important, 
pions are shown separately in Fig. [l](b). Thus, a pion chemical potential of nearly 80MeV 



may be acquired. Similar values have been found previously |T5| , |24fl . See 0] f° r a discussion 
of how inelastic reactions can reduce these values. 

Now we return to the EOS. The energy density as a function of temperature, with and 
without the chemical potentials, is shown in Fig. [| The principal observation is clear: with- 
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FIG. 2: The energy density as a function of temperature with and without chemical freezeout for 
the SPS (s/n B =42). 

out chemical freezeout the energy density drops very rapidly as a function of temperature, 
since the total number of particles drops. For later discussion, the temperature is tabulated 
as a function of energy density with and without chemical potentials, in Table |. However, 
to find the hydrodynamic solution the most important quantity is not the relationship be- 
tween energy density and temperature, but the relation between the energy density and 
pressure - the EOS. Fig. || shows the pressure, the squared sound speed c 2 s1 and the entropy 
density versus the energy density along the adiabatic path. Along the adiabatic path, these 
quantities are all related because 

/ dp\ 



^ ' {ni/s} 

ds\ 

dels i \ 



p + e 



(The second relation follows by noting that ( ^ ) = ( 

V " e / m /s V 



ds 

drii / e 



ii ( d±\ 



(12) 
(13) 

, solving for 

(^H and using thermodynamic identities). Unlike energy density versus temperature, 
pressure is scarcely modified by the chemical potentials. Differentiating the pressure, we see 
that c 2 s is reduced somewhat at low temperatures when chemical freezeout is incorporated. 
Integrating to find the entropy, we see that the entropy as a function of energy is nearly 
identical with and without chemical freezeout. 
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Energy Density (GeV/fm 3 ) 


T(MeV) with fi n > 


T(MeV) with = 


0.364 


160 


163 


0.249 


140 


153 


0.158 


120 


142 


0.122 


110 


135 


0.091 


100 


129 


0.067 


90 


122 


0.047 


80 


115 



TABLE I: Temperature versus energy density, with and without chemical freezeout at the SPS 
(s/n B = 42). 

From the point of view of dynamics, this means that the hydrodynamic solutions, with 
and without chemical freezeout, are nearly the same when expressed as a function of energy 
density. The flow velocities on a freezeout surface of constant energy density are independent 
of whether or not chemical potentials are included. However, the temperature on that 
freezeout surface depends dramatically on chemical freezeout, as can be seen from Table |. 
To make this point clear, the hydrodynamic solutions are illustrated in Fig. |] with and 
without chemical potentials. The solutions are similar; the small differences can be traced 
to small differences in the speed of sound. 

Often in hydrodynamic simulations not incorporating chemical freezeout, the freezeout 



surface T « 130 MeV is taken [17], [18| . Using Table 0, this roughly corresponds to a temper- 
ature of T « 100 MeV. Although this temperature seems low, it is not out of keeping with 
the phenomenological freezeout temperatures extracted from thermal fits to radial and 



recent elliptic flow data [2?J. The extent to which this low freezeout temperature is seen in 
a hadronic cascade is addressed in the next section. 



III. RQMD AND CHEMICAL FREEZEOUT 

Returning to hydro+cascade, we can study the extent to which chemically hydrodynamics 
reproduces the dynamic response of the cascade. An important observation is that with the 
assumption of chemical equilibrium, the hydrodynamics is independent of the dominant 
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FIG. 3: The (a) pressure, (b) sound speed squared and (c) entropy density as functions of the 
energy density with and without chemical freezeout at the SPS {s/ub =42). The analogous curves 
at RHIC are only slightly different. 

reactions in the cascade. With chemical freezeout, reactions such as ttN — > A — > ttN, 
are encoded into the conservation laws and ultimately into the chemical potentials which 
affect the spectra. Baryons are accelerated only by reducing the kinetic energy of the pions. 
In a hydrodynamic language, the flow velocity is increased only by rapidly decreasing the 
temperature. 

One of the problems with the hydro+cascade approach JET]] is the sensitivity to the 
switching temperature. To make a smooth transition from hydrodynamics into the cascade, 
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FIG. 4: The hydro dynamic solution (a) with and (b) without chemical freezeout at the SPS 
(PbPb, y / s=17 GeV A, b=0 fm, s/n^=42). The thick arcs show contours of constant energy 
density. The first contour indicates the start of the mixed phase, eg. The next three contours 
indicate energy densities corresponding to temperatures (a) T = 160, 120, 80MeV with chemical 
freezeout and (b) T=163, 142 , 115 MeV without chemical freezeout, see Table |. The thin lines 
shows contours of constant transverse fluid rapidity, y^ = 0.1, 0.2, . . . , 0.7 . Walking along the thick 
arcs, the arc is divided into solid and dotted segments. 20% of the total entropy passing through 
the arc passes through each segment, (yx) denotes the mean transverse rapidity (weighted with 
entropy) on the arc. 

the same conservation laws have to be implemented in both approaches. In Fig. |, we adjust 
the switching temperature from T switch = 160 MeV to T switch = 117 MeV , and compare 
spectra in which chemical freezeout is incorporated, fi n > 0, to spectra in which chemical 
freezeout is not incorporated, fj, n = 0. Ideally, the spectra should be insensitive to T switch . 
First, notice that if chemical freezeout is not incorporated the yields of K and p are re- 

_Mk__m jl m n m n 

duced by factors of e T T = and e T T = , or 3 and 7 for temperatures of T th 120 MeV. 
Furthermore, the final flow is too strong, since both the the mass energy and the kinetic 
energy are converted into flow. Second, notice that once chemical freezeout is incorporated, 
the sensitivity to the switching temperature is small. This is because the freezeout energy 
density is higher; therefore it is "as if" the freezeout temperature were higher. Close inspec- 
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FIG. 5: Sensitivity of particle spectra at the SPS (PbPb, yfs=Y7 GeV A, b=6 fm, s/ne=42) to 
Tswitch with chemical freezeout, p, w > and without, fi n = 0. The spectra are for ir~ , p, and p. 
The filled histograms are T sw i tc h = 160. The other curves are for T sw i tc h = 117. 

tion of the pion spectrum shows that the final spectrum with chemical freezeout is too cool 
relative to the Hydro+cascade model. This rapid cooling will be addressed below. 

A naive explanation for the insensitivity to T SW i tc h is that the decrease in temperature 
is compensated by flow. Fig. || studies this explanation in detail. Hadronic rescattering 
increases the mean M? of the nucleons while slightly decreasing the mean My of the pions. 
Although these general trends are reproduced by the hydrodynamics, with the hydrodynam- 
ics the pion cooling is too rapid and the nucleon acceleration is slightly too large. These 
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FIG. 6: The effect of hadronic scattering on the proton spectrum at the SPS (PbPb, y/s=17 GeV 
A, b=6 fm, s/n#=42) with and without chemical freezeout. The temperature T sw it c h, where the 
switch from hydro to cascade made is given in the figure. H(ydro)+RQMD denotes the spectra 
with subsequent hadronic rescattering in the cascade. H(ydro) Only denotes a spectrum in which 
only resonance decays but no rescattering have been accounted for. 

facts may be gleaned from a very close inspection of the four curves in Fig. || However, 
the differences are all small and it is difficult to tell the difference between chemically frozen 
hydrodynamics and free streaming from spectra alone. 

Instead of spectra, we study the elliptic flow when chemical potentials are included. Fig. [7] 
shows the sensitivity of elliptic flow to T switc h. The elliptic flow remains fairly sensitive to 
the switching surface and increases by 25% as the switching temperature is lowered. The 
hydrodynamic contribution increases while the RQMD contribution decreases as T switc h is 
lowered. The fact that the Hydro+RQMD curve is flatter than the Hydro Only curve 
indicates that hydro+cascade model is partially successful. 

Additional information is gained with Fig. [8|, which examines the elliptic flow spectrum, 
i> 2 (pt)- The rapid cooling of the hydrodynamic response with chemical freezeout (T=117 
fx n > 0, Hydro+RQMD) is seen as an increase in elliptic flow at modest px, relative to the 
normal curves (T=160 Hydro+RQMD). By comparing the T = 160 MeV Hydro Only and 
Hydro+RQMD points, we see that in the cascade the pions cool only slightly. Nevertheless, 
the general trend is the same for both the hydro with chemical freezeout and the cascade 
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FIG. 7: V2 of pions at the SPS (PbPb, y/s=17 GeV A, b=6 fm, s/ns=42) as a function of T sw i tc h 
with chemical freezeout, fi n > 0, and without, ^ = 
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FIG. 8: The elliptic flow spectrum of pions at the SPS (PbPb, ^/s=17 GeV A, b=6 fm, s/n B =42) 
as a function of T sw i tc h with chemical freezeout fi n > 0, and without, /x^ = 0. 
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with chemical freezeout. 

We now turn to the distribution of pions at freezeout. The density of pions at freezeout 
is approximately constant as a function of impact parameter and beam energy |TJJ. Now 



we interpret this density in thermal terms. In the cascade, the density of pions at freezeout 
is the same as the density of pions in a hadron gas with T rs 110 MeV and ^ ~ 70MeV. 
Notice Fig. [| which shows the mean emission time as the switching temperature is lowered. 
As with the spectra, the mean emission time after cascading is insensitive to the switching 
temperature. The system of hadrons introduced into the cascade expands until the density 
reaches a certain value and subsequently breaks up. Since the number of hadrons is con- 
served during the hydrodynamic evolution, the final breakup does not depend on where the 
hydrodynamics is stopped and where the cascade begins. 

We now try to measure the RQMD pion density at freezeout. The cascade does not have 
a freezeout surface, as is normally assumed in an idealized hydro picture. Rather, particles 
are emitted per unit time and volume, i.e. , ~ , , . We define the freezeout density in the 

* ' rdrj ax ay dr J 

transverse plane as 

n RQMD = 1 f dxdydT dN 

-kR^ J^/x 2 +y 2 <R r arj ax ay dr 

where R Q is taken to be 3fm. This density, n^ MD , with and without hadronic rescattering 



is shown in Fig. [LC]. Of course, without rescattering, the density of pions should reflect the 
density of pions on the freezeout surface. With the EOS used in this work, this number 
density is simply the zero-th component of j£ or 

Jh* — n K + 2n p + n A (15) 



This thermal density of pions is also shown in Fig. |10| and gives a reasonable description of the 
Hydro Only points. When rescattering is included, the density decreases until a constant 
value ~ 0.1 fm 3 , is reached. Comparing the thermal model curves to the Hydro+RQMD 
points, we see that the freezeout densities are equal for T « 100 — 110 MeV and /x w ~ 80 MeV. 
Although this temperature is low, it is not out of keeping with the numbers extracted from 



hydrodynamic fits to the data [2§, ^ . 
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FIG. 9: The mean emission time < r > at the SPS (PbPb, y/s=Yl GeV A, b=6fm, s/n B =42) as 
a function of T sw n c h with chemical freezeout (j, n > 0, and without chemical freezeout, fi n = 0. 




FIG. 10: The pion density at freezeout in RQMD (see Eq. [l4|) as a function of T sw i tc h at the SPS 
(PbPb, y / s=17 GeV A, b=6 fm, s/n j B=42). The Thermal Model curves indicate the density of 
pion number (see Eq. |l5|) in a hadron gas, i.e. n n + 2n p + ... . 
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IV. SUMMARY OF CHEMICAL FREEZEOUT 



We have constructed an EOS consistent with chemical freezeout. The relationship be- 
tween pressure and energy density is approximately the same as the standard EOS. However, 
these two EOS exhibit dramatically different relationships between energy density and tem- 
perature. At a practical level, this means that the results from a hydrodynamic calculation 
without chemical freezeout may be consistently converted to a calculation with chemical 
freezeout by simultaneously adjusting the yields and the temperature to keep the energy 
density constant. Table | gives the required conversion factors. 

We used this chemically consistent EOS to study the hydrodynamic + hadronic cascade 
model, H2H |[L4|| . Since, chemical freezeout is incorporated via the hadronic cross sections 
into the cascade, the same conservation laws are implemented in both the hydrodynamics 
and the cascade. With this congruence, we examined the model sensitivity to the switching 
temperature from hydro to cascade, T switc h. The spectra are insensitive to the switching 
temperature for 120 MeV < T SW i tc h < 160 MeV. On the other hand, elliptic flow at the 
SPS remains mildly sensitive to T switc h even when chemical freezeout is incorporated. The 
insensitivity of the model's results to T switch partially validates the hydro + cascade approach. 

An important feature of a chemically frozen EOS is that the principal reactions in the 
cascade, e.g. ttN — > A — > ttN, are encoded into the hydrodynamics through the chemical 
potentials. Thus, in the cascade, the nucleons are accelerated only by reducing the kinetic 
energy of the pions. With a chemically frozen EOS, this qualitative feature is reproduced by 
the hydrodynamics. However, as might be expected, the cooling and attendant acceleration 
are larger in the hydrodynamics than in the cascade. When the cascade is replaced with the 
hydrodynamics, the pion spectrum is slightly steeper while the nucleon spectrum is slightly 
flatter. Pion cooling may also be seen with ^(Pt)- Here again, the small changes in ^(pt) 
due to hadronic rescattering are strongly magnified by the hydrodynamics with chemical 
freezeout. 

By incorporating chemical freezeout, many of the qualitative features of the cascade's 
evolution are reproduced with the hydrodynamics (see ]30|] for further investigations of this 
point). Thus, a thermal interpretation of the cascade's response is partially justified. Within 
RQMD, we find the density of pions at freezeout to be n n « 0.12 fm~ 3 , which corresponds 
to T « 110 MeV and /x w ~ 70 MeV in thermal terms. Precisely these parameters have been 
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extracted from thermal fits to hadronic data. Therefore, an understanding of the properties 
of a chemically frozen hadronic gas helps to bridge the cascade and thermal descriptions of 
the final stages of the heavy ion collision. 
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